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Abstract 



. Several nuclear reactions involving the ^Be nucleus, not included into the 

^ i standard model of the pp-chain, are discussed. A qualitative analysis of their 

' possible influence on the fate of the ^Be in solar plasma and of their role in 



the interpretation of the solar neutrino experiments is given. As an example, 
' the reaction rate of the nonradiative production of ^B in the triple collision 

p + e~ + ''Be ^B+e~ is estimated in the framework of the adiabatic 
approximation. For the solar interior conditions the triple collision reaction 
^ ^ rate is approximately 10~^ of that for the binary process p + ^Be — > ^B + 7. 

§ ; [PACS numbers: 21.45. +v, 95.30.-k, 97.10.Cv] 

> : 

• 1— I , 

^ ■ I. INTRODUCTION 



The so-called standard model of the sun has been developed during several decades by 
a collective effort of many researchers. At present it is an elaborate theory which accounts 
for many, though not all, observable characteristics of the sun. One of the exceptions is the 
flux of neutrinos which is several times higher than what is actually seen in experiments 
(for details on the neutrino problem see, for example, Ref. [0). There are several possible 
reasons for this discrepancy which are associated either with the nuclear processes involved 
or with the neutrino propagation towards the earth. 

The solar neutrinos emerge from nuclear reactions which occur in plasma consisting of 
light nuclei (mainly protons) and electrons. The sequence of these reactions begins with the 
proton-proton collision and is, therefore, called the pp-chain. In the standard model it is 
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assumed that this chain consists of the following reactions 



p + p ^H + e+ + z/ 

p + p + e" ^H + z/ 

p + ^H 2He + 7 

^He + ^He ^ ^He + p + p 

3He + ^He ^ ^Be + 7 

and then the ^Be is destroyed in two ways 

e" + ^Be ^ '^Li + z/ p + ^Be ^ ^B + 7 

p + ■^Li ^Be + 7 *^B ^ ^Be* + e+ + 

^Be ^ ^He + ^He ^Be* ^He + ^He . 



The fate of ^Be in this chain is of special interest since a combined analysis of all experiments, 
in which the neutrino flux from the sun was measured, has come to a paradoxical conclusion: 
the production of ''Be nuclei (more precisely, the flux of neutrinos due to ''Be reactions) must 
be strongly suppressed or even be negative P,^. This implies that something is wrong either 
in the standard model or in the experimental data. 

In this work we present a qualitative as well a quantitative analysis of the nuclear reac- 
tions involving ''Be, which may contribute to a more adequate description of the behavior 
of this nucleus in the sun. Our approach differs from the standard model in that we take 
into account also the triple collisions in solar plasma. In contrast the standard pp-chain 
consists of a sequence of reactions of two-body initial states, the only exception being the 
ppe~ collision. 

There are two essential differences between the nuclear reactions caused by binary and 
triple collisions. They can be classified as kinematical and dynamical. The former are 
related to the selection rules prevailing in the two- and three-body reactions while the 
latter stem from the inter-dependence of different nuclear processes. Thus, some binary 
nuclear reactions suppressed by conservation laws (concerning angular momentum, parity, 
isospin etc.) can be eased due to the presence of a third particle. Therefore, the three-body 
mechanism of such reactions which is less restricted kinematically, may play a significant role 
in the nuclear burning in stellar plasma where the probability of triple collisions could be 
quite high due to the high density of matter. The dynamics of the three particle motion may 
also lead to a completely different physical picture: Processes considered as independent in 
binary collisions become dependent when triple collisions are considered. For example, the 
e' + ^Be and p + ^Be processes, become dependent to each other when the e~ + p + ^Be 
collision is taken into account. 

The paper is organized as follows. In Sect. II we discuss the role of the three-body 
collisions qualitatively. In Sect. Ill we describe our formalism and in IV and V we outline 
the procedure employed to evaluate the various ingredients used to obtain the reaction rate. 
In Sect. VI we present our results and conclusions. Some details concerning the derivation 
of the reaction rate formula are given in the Appendix. 
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II. BINARY AND TRIPLE COLLISIONS 



According to the standard model of the sun the ''Be nucleus, produced in the pp-cycle, 
is destructed via the following two binary reactions 

e- + ^Be^ ^Li + z/, (1) 

p + ^Be^ 8B + 7. (2) 

Due to the high density of the plasma, proton and electron can always be found in the 
vicinity of ''Be and form a three-body initial state 

p + e" + ^Be) , (3) 



which can give rise not only to the above binary transitions (|l|) and (0) but also to the 
processes 

r ^Li + p + 1/ 

p + e" + ' Be — ^ < *B + 7 + e" . (4) 
i SB + e- 

The first reaction corresponds to the weak transition (|T]) provided the nearby proton does 
not participate in the reaction. Similarly, the second of these reactions corresponds to the 
radiative capture (||) if the electron is a spectator particle. It is obvious that the three 
processes are not independent, since they are generated from the same initial state. 
Therefore, the neutrino fluxes related to ^B and ^Be nuclei are not independent either. This 
means that the procedure of comparison of the Homestake and Kamiokande experiments, 
which gave patently wrong negative amount of ^Be, could be incomplete and should be 
reconsidered. 

The question then arises on how strong is the mutual dependence of the processes (||). 
This question can be formulated in an alternative way: To what extent the two-body sub- 
systems, involved in (^), can be considered as independent? Formally such an independence 
would imply that the total wave function of the initial state @ could be written as a direct 
product of the wave functions of the e~+^Be, p+''Be, and p+e~ subsystems, 

p + e" + '^Be) ^ e" Be) (g) p Be) (g) p + e~) . (5) 



The possibility of such a factorization for three charged particles was investigated in Refs. 
0J^, where it was shown that the total wave function reduces to a product of the type 
only when all three particles are far away from each other. Such a configuration, however, 
can not contribute much to the transitions for they are caused by the strong and 
weak interactions which vanish at large distances. The description of the processes as 
independent needs therefore an additional substantiation. 

There are also other processes which, in principle, can change the balance between the 
neutrino fluxes from ^Be and ^B without changing the abundance of ^Be nuclei. Two such 
examples are 

e" + e" + '^Be ^Li + e" + z/ , 
p + p + ^Be ^B + p . 
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In this article, however, we investigate only the processes leading to the formation as 
described by (1). This requires the reaction rate of the process 

p + e" + ^Be -> *B + e" (6) 

which is discussed next. 



III. REACTION RATE 

To estimate the reaction rate of (|^) we assumed that the nucleus ^Be can be considered 
as a paricle and, therefore, the ^B nucleus is a bound state of ''Be with the proton. Of course 
such an approximation is inadequate if a full description of the spectrum and the structure 
of *B is required. Since, however, its ground state is only 0.138 MeV below the (p+'^Be)- 
threshold (see Fig. 1) this cluster representation should describe it fairly well. In the 
problem under consideration the other states of ®B are not needed since the collision occurs at 
ultra-low energies of the order of ~1 keV (at temperatures of ~ 10'' °K) and thus transitions 
via excited states can be safely neglected and the formation of ®B proceeds directly from the 
continuum (p+^Be) through a kind of Auger transition in which the electron carries away 
the excess energy. 

The Hamiltonian, which governs the evolution of the three-body initial state (^, can be 
written as 

H = H^ + hQ + Vn + Ve, (7) 

where Hq and are the kinetic energy operators associated with the Jacobi variables r and p 
respectively (see Fig. 2), is the p^Be potential, which includes strong and Coulomb parts, 
and Ve is the sum of the e~p and e~''Be Coulomb potentials. The amplitude, describing 
various transitions in this three-body system, is proportional to the T-operator obeying the 
Lippmann-Schwinger equation 

T{z) = K + K + (V; + K) ^ -T{z) . (8) 

The initial asymptotic state (^ is defined by the p^Be relative momentum p, the total 
nuclear spin s and its third component m^, and by the electron momentum with respect to 
center of mass of the nuclei k. The final asymptotic state consist of the bound state of a 
proton and ^Be, with spin a = 2 and third component m^, and a freely moving electron with 
momentum k' which differs from k not only in its direction but also in its absolute value, since 
the electron picked up the energy released in the nuclear process. We note that we neglect 
the spin of the electron since we assume that its interaction can be satisfactorily described by 
purely coulombic potentials without spin-orbit or other forces involving magnetic moments. 
Therefore, the reaction (^ can be described as a transition between these asymptotic states, 

KZf,k') . (9) 

The corresponding three-body reaction rate 71^ is given by (see Ref. 0) 



p, s. 



TUs'-, k 
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7^3(p, s, rris, k ^ am,, k') = -^^Ej - E^) \{^lZ'l\ k'\ T |p, s, m,; k)'\ n.n^n,^, , (10) 

where the (5-function secures the energy conservation while Ueripm^f. is the product of the 
densities of electrons, protons, and ^Be nuclei in the plasma. The coefficient I/Att"^ in Eq. 
(p!0|) stems from the following normalizations of the wave functions 



^bound 
ania 



k 



k') = {2Tif5{k-k') 



Both possible values of the p^Be spin, s = 1 and s = 2, as well as all its orientations are 
represented by different p''Be-pairs in the plasma with equal probabilities. The momenta p 
and fc, however, are distributed at the temperature 9 according to Maxwell's law 



Nj:{e) = {2nfiKB0y^ exp (^-k^/2fiKB0 



where kb is the Boltzmann constant and the reduced masses m and fi are defined by m ^ = 
+ m^gg and = + {nip + rm-Q^)'^. 
We need the average reaction rate due to all possible initial states with different p, s, m^, 
k and the total rate of transition to the final states with any allowed ma and k'. This can 
be obtained by averaging the rate ([T0|) over all initial and summing it over all final quantum 
numbers. The resulting reaction rate depends only on the temperature of the plasma 

<n,{e)>= Yl fd'k' fd^p fd'kNpie)Nj;{e)ns{p,s,m,,k^ama,k') . (11) 

In what follows we consider the various approximations and ingredients needed to evaluate 
the reaction rate by means of this general formula. 



IV. ADIABATIC APPROXIMATION 

The average kinetic energy of the particles in the plasma, < E^^"' >~ ^bS ~ 1 keV, is the 
same for nuclei and electrons, but the velocity of an electron is three orders of magnitude 
higher than that of a proton or ''Be. Therefore, while the proton approaches ''Be very slowly, 
the electron dashes nearby picking up the energy and leaving the heavy particles in a bound 
state. Due to this we can significantly simplify the problem by treating the relative p^Be 
motion adiabatically. Such an approximation means that the electron bypasses the heavy 
particles very fast and thus it has no time to 'observe' any changes in their space position. 
Moreover, to participate in such three-body reaction, the heavy particles must be close to 
each other when the electron arrives, because they are too slow as compared to the electron. 
Hence, while the electron starts from its asymptotic state the heavy particles have 
already interacted to each other via the potential Vn and formed the two-body scattering 
state 
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\v,s,m,) ^ ^^S^l) . (12) 

Therefore, instead of the transition (|^) caused by both Vn and Vg, in the adiabatic approxi- 
mation we may consider the transition 



■^scat 1 
psnis ' 



Kz:^k') , (13) 



where the interaction Vn is taken into account by (|T^). 

Formally, this leads to some simplifications in the Eq. (^. Indeed, since v|/**c<^* and v]/'"'""'^ 
can be considered as 'excited' and ground states of the nuclear target respectively, we can 
describe the transition (|T3|) using another T-matrix defined as 

= K + V ^ f{z) , (14) 
z — ho — Ha 

where = i^o + ^ is the total Hamiltonian of the nuclear subsystem. Within the above 
approximation we have 



■^Tjbound j^/ 



T 



p, m.; k) ^ (^iZ:\ k' f ^^.t,, fc> . (15) 



It is known (see Ref. [§]) that the adiabatic approximation is equivalent to the closure 
approximation that is when all excited states of the target are assumed to be degenerated 
and the Hamiltonian is replaced by a constant, viz. 

Ha = J2^n\ipn >< i^n\ + f dE E^je >< V'sl 

^^ojl^l^n >< ^nl + j dE\^E><'4'E\^ =^0 • 

Therefore, in the adiabatic approximation T ^ T°, where is defined by the following 
equation 

f\z) = v; + v; — ^^T\z) . (16) 

Z — tlQ — Cq 

In our case Sq is the ground state energy of the nucleus, Sq = — 0.138 MeV. Substituting 
into Eq. (|16D the total energy z = E^^^ + Sq, where E^^^ is the kinetic energy of the electron 
with respect to the center of mass of the heavy particles, we obtain the Lippmann-Schwinger 
equation describing the scattering of an electron from fixed centres. For the solar plasma 
electrons e'^/hv < 1, which is a sufficient condition to treat Coulomb interactions in Eq. (|TB[) 
perturbatively, i.e.. 

Furthermore, the average potential energy of the electron is of atomic order of magnitude, 

< Ve >~ 10 eV, while its kinetic energy in solar plasma is two orders of magnitude higher, 

< E^^"" >~ 1 keV, which implies that the above iterations should converge very fast. There- 
fore, we may retain only the first (Born) term 0. 
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Finally, the exact formula (|TOD for the reaction rate can be replaced by the following 
approximate one 



arria 



k') 



47r2 



5{Ef - E,) 



bound 



k' 



^scat 1 
psrris 5 ^ 



(17) 



which can be considered as a three-body generalization of the distorted wave Born approxi- 
mation (DWBA) which is widely used in the theory of the two-body nuclear reactions. The 
explicit formula for the reaction rate, which involves (^) ^CTm"'^(^ described in the 
next section, is derived in the Appendix. 



V. NUCLEAR WAVE FUNCTIONS 
A. Scattering state 

As mention earlier at the ultra-low energies considered (~lkeV) the collisions between 
protons and ''Be nuclei are dominated by the S-wave state. Therefore, in the partial wave 
expansion of the corresponding scattering wave function we may retain only the term with 
zero angular momentum 

"^fsti^ = - E < ^/^i|/^2km, > xUlT , (18) 

where and x^^^ are the spin functions of the proton and ^Be nucleus. The S'-wave radial 
function ^^^^(r) is a regular solution of the radial equation 

1^ + - 2m [Vc{r) + Vs{r)] | <'^*(r) = (19) 

with the asymptotic boundary condition 

uf^^Ur) — y sin \pr — ri\n (2pr) + Sq] , 

where t] is the Sommerfeld parameter. The potentials Vc and Vs in Eq. (p!9| ) describe 
the Coulomb and strong p-^Be interactions respectively. In order to take into account 
the nonzero size of ^Be nucleus we assume that it is a uniformly charged sphere of radius 
R = 2.95 fm. Thus the potential Vc is 

f 1/r , r>R 

Vc{r) = 3 (20) 

[2^-2^' 

The strong potential is assumed to have a Woods-Saxon form [|10| 



^'^''^ ^ ~l + exp[(r-i?)/ajv] ^^^^ 



with Vq = 34MeV, = 0.52 fm, and the same nuclear radius R. 
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B. Bound state 

The ground state of ''Be has the quantum numbers 3/2^. Therefore to construct the 
ground state of the Boron isotope ^B with angular momentum 2 and positive parity, in the 
two-body p^Be-model, the angular momentum I of p^Be relative motion must he I = 1. 
The wave function of ^B can then be written as 

KZfi^ = - E < lf^4f^2\sm, X sm,lm\2m^ > x^X^f • (22) 

The radial wave function 0''o™'^(^r) is a square-integrable solution of the equation 

^ + 2m^ - ^ - 2m [^^(r) + V^(r)] | 0"(r) = , (23) 

where the potential V'^ has the same form as ( pT]) but the strength parameter is slightly 
modified, = 32.62 MeV to reproduce the bound state of *^B, So = -0.138 MeV. 



VI. RESULTS AND DISCUSSIONS 

The average production of ®B nuclei per unit volume per second via the triple collisions 
(§) at a temperature 6 can be written as (in units cm~^ sec~^) 

< 7^3(^) >=< S(0) > nenpnTBe 

where < S(6') > is analogous to the so-called transport cross-section < > used in the 
theory of two-body reactions [|11[]. Both < T,{6) > and < > are the reaction rates 



corresponding to unit densities of particles. The explicit form of < S(^) > is given in the 
Appendix. 

The average production of ^B via the radiative capture (|^) is (in units cm^'^ sec^^) 

< '^2(6') > = < (T2V > npWTBe 

Therefore, the ratio of the triple to binary rate depends only on the density of electrons and 
the temperature 

< ^3(g) > _ <m> ^ .24) 

We calculated this ratio using Ue = lOOA^yi cm~'^ where Na = 6.022 x 10^^ is the Avogadro 
constant. This density corresponds to solar plasma conditions [0. The results of our calcu- 
lations for several values of 6 are given in Table 1. For comparison, this table contains also 
the absolute values of the reaction rates of the radiative and nonradiative processes (0) and 
(I) calculated for = mBe = Na- The data for the radiative proton capture by ^Be are 
taken from Ref. It is seen that the triple collision (^ plays, apparently, a minor role 
similarly to the nonradiative capture reaction e~ + p + ^ '^He + e~ |jl2[|. However, at the 
early stages of the universe when ne/NA ^ 100, this reaction, might have been significant. 
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The temperature dependence of these reaction rates, normahzed to the same value at ^ = 
14 X 10^ °K, are shown in Fig. 3. It is seen that the curves for the radiative and nonradiative 
processes practically coincide. This is not surprising, for in both cases the temperature 
dependence is determined mainly by the same exponential function stemming from the 
Maxwellian distribution of the p''Be momentum p (see formula (|A9| ) in the Appendix). 
Indeed, a change of the temperature shifts the maximum of the distribution which changes 
the most probable energy of p^Be collision. This in turn affects the values of '^^smS'^ near 
the point r = 0. Since the Coulomb interaction between the proton and ^Be is repulsive 
and the collision energy is very low, the scattering wave function at small distances depends 
on the temperature very strongly. This dependence manifests itself via the corresponding 
behaviour of the reaction rates because both the radiative and nonradiative transitions occur 
at short distances and their rates depend on the overlapping of ^E'|§ms(^) with the bound 
state wave function \l/^'^^'^(r). 

Of course, the energy of the electron also depends on the temperature. This however has 
not much influence on the reaction rate because the electron is attracted to the area of the 
reaction and its wave function is always large at small p. An additional increase or decrease 
of it makes, therefore, no much difference in contrast to the repulsive interaction. 

This reasoning leads us to yet another conclusion. A reaction with a three-body ini- 
tial state and with all particles positively charged, must have very strong temperature de- 
pendence. Thus one should expect quite different 6'-dependence for the nonradiative ^B 
production with two protons in the initial state 

p + p + ^Be^ ^B + p. 

Since in this case the adiabatic approximations is not expected to work well, a more rigorous 
treatment of this reaction, for example, using the Faddeev formalism, is required. 
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APPENDIX A: EVALUATION OF THE TRANSITION RATE 



The potential Vg describing Coulomb interaction of the electron with the proton and ^Be 
in the coordinate space is 

Vg{p,f^ = - (Al) 

\p + jjr\ \p — ar\ 

where 

a = , p = 1 — a . 

m7Be + rup 

The Jacobi vectors p and r* are depicted in Fig. 2. According to (plTf ), the matrix element 
between the initial and final 'distorted' waves provides the amplitude of the transition (p!3D 



M{p, s, m„ k am^, k') = (vj/^— ^ k'\ K k) = J df^lZ:'*{miq, ^n^li^ ' 

(A2) 

with 

^ 47r ^ ^ 

Ue{q,r) = [4exp(— ■ r) + exp(iaq- r)] , (A3) 

where q = k — k' is the momentum transferred to the electron, and Ue{q,r) is the Fourier 
transform of the two-center potential (|A1| ) over the electron variable p. Performing the 
angular integrations in (|A^), we obtain 



z(47r)^ 

PQ m 

(A4) 

where the radial integrals are given by 



oo 



Ril,q,p)=l 0"(r)<-*(r)ji(7gr)dr (A5) 







with ji being the Riccati-Bessel function As can be seen from Eq. (|A4|) and the 

definition of q all partial waves of the electron motion contribute to the transition amplitude 
M{p, s, rris, k — > am^, k'). 

By averaging over the initial and summing over the final spin orientations we obtain 

^ \M{p,s,ms,k^am,,k')'' = ■^^^\R{a,q,p)-AR{l3,q,p)f . (A6) 

Then the reaction rate (ITTF) can be written as 



< 7^3(^) >=< S(^) > rienprnee , (A7) 

where 
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J2 



(A8) 



\R{a,q,p)-AR{l3,q,p)f . 



The integration over k' in the above formula can be performed using q=k' — k instead of k' 
in the integration. Then the 5-function facihtates the integration over the sohd angle of q. 

Finally, the unit density reaction rate < S(6') > is expressed in the form of a two- 
dimensional integral 



< m > = E 



10(27r)4^ 
2s + 1 ^2 7o 



dpexp 



X £ -exp 



roo 

7^ 



2eq^ I 2 2ix ° 



(A9) 



[Q{a,p,q)-AQ{(5,p,q)f 



with 



Q{a,p, q) = (j)'"''''"^{r)jx{aqr)ul^"'\r) exp f ^ j dr . 
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TABLES 





< > Na 
( cm^ mol^-"^ sec^-*^) |Tl]l 


< > Nl 
( cm'^ mol^i sec^i) 


< n^ie) > / < 7^2(^) > 


14 


1.73 X 10-12 


1.83 xlO-i*^ 


1.06x10"^ 


15 


4.34 X 10^12 


4.43 X 10^16 


1.02x10-^ 


16 


10.1 X 10^12 


9.80 xlO^i*^ 


0.97x10"^ 


18 


44.7 xl0'i2 


42.0 xlO-i*^ 


0.94x10"^ 


20 


161.0x10^12 


143.3x10^1^ 


0.89x10"^ 



TABLE I. Temperature dependence of the reaction rates for the radiative and nonradiative 
formation of '^B in solar plasma. 
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FIGURES 
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FIG. 1. Spectrum of nucleus. Energy levels are in MeV. 
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FIG. 2. The Jacobi vectors. 
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FIG. 3. Temperature dependence of the binary and triple reaction rates given in Table 1 and 
normalized to the same value at ^ = 14 x 106 °K. 
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